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0^ , Abstract. We apply the method of nonlinear steepest descent to compute the long- 

' time asymptotics of the Camassa-Holm equation for decaying initial data, completing 

I previous results by A. Boutet de Monvel and D. Shepelsky. 

^ : 

! 1. Introduction 

. 

■ In this paper we want to study the long-time asymptotics of the Camassa-Holm 
(CH) equation, also known as the dispersive shallow water equation, 

I— I ■ 

C/!5 I Ut + 2>iUx-Utxx + 3uUx = '^UxUxx+UUxxx, t > 0, XGR, (1.1) 

I where u = u{x,t) is the fluid velocity in the x direction, x > is a constant related 

■ to the critical shallow water wave speed, and subscripts denote partial derivatives. It 
I— I, was first introduced by R. Camassa and D. Holm in [7J and R. Camassa et al. [8j as a 

model for shallow water waves, but already appeared earlier in a list by B. Fuchssteiner 
^ ■ and A. Fokas [19]. More on the hydrodynamical relevance of this model can be found 

^-H I in the recent articles by R. Johnson [22] and A. Constantin and D. Lannes [12] . With 

Cn ■ w := u - Uxx + (1-2) 

o ■ 

' called the "momentum", equation (II. ID can be expressed as the condition of compati- 

Q ■ bility between 

S: i(-/" + i/)=A/. (1.3) 



(N 



that is, 

dtdxxf = dxxdtf 

is the same as saying that (jl.ip holds. Equation (jl.3p is the spectral problem associated 
to (jl.ip . In particular, the CH equation is completely integrable and can be solved 
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via the inverse scattering method. Correspondingly, we consider real-valued classical 
solutions u{x,t) of the CH equation (jl.ip . which decay rapidly, that is, 

m» /■(l + |.|)'«(W.,()-.| + |„.(.,<)l + K,(..OI)d.<co (L5) 

0<t<T Jjj 

for all T > and some integer / > 1. Moreover, we will assume 

w{x,0)>0 (1.6) 

throughout this paper. Then u exists for all times t > with w{x,t) > (for exis- 
tence of solutions we refer to A. Constantin and J. Escher [T^ and the discussion in 
A. Constantin and J. Lenells [E], see also [9]). 

The aim of this paper is to establish the long-time asymptotics of such solutions 
using the nonlinear steepest descent method from P. Deift and X. Zhou [16] which was 
inspired by earlier work of S. Manakov [27] and A. Its 21]. More on this method and 
its history can be found in the survey by P. Deift, A. Its, and X. Zhou [17j. 

The starting point for this method is the representation of a solution of the nonlinear 
equation under consideration in terms of a solution of an associated Riemann-Hilbert 
problem. 

Recently, A. Boutet de Monvel and D. Shepelsky have used the inverse scattering 
approach to the CH equation, based on the construction and analysis of an associated 
matrix Riemann-Hilbert problem [Sid]. The analysis, via the nonlinear steepest descent 
method, of a vector oscillatory RH problem derived from the original matrix- valued 
problem, allowed [Sj distinguishing four main regions in the (x, i)-half-plane, where the 
leading asymptotic terms were qualitatively different: a solitonic sector, two sectors of 
(slowly decaying) modulated oscillations, and a sector of rapid decay. 

Here we want to simplify the original approach, deriving the vector Riemann-Hilbert 
problem directly from scattering theory for the underlying Sturm-Liouville operator. 
We notice that the matrix and vector RH problems, being closely related, have specific 
features concerning, particularly, the uniqueness issue. For the matrix problem, we 
refer to [5], whereas the uniqueness for the vector problem is addressed in detail in the 
present paper, see Section [3] below. 

At the same time we want to provide complete proofs including the effects of solitons 
in the sectors of decaying oscillations (formulas given in [5j for the oscillatory regions 
are true in the solitonless case only) and the error estimates in terms of decay of the 
initial condition. Moreover, we will make use of some simplifications to the nonlinear 
steepest descent method recently given in H. Kriiger and G. Teschl [23] respectively 
K. Grunert and G. Teschl [20] . 

The asymptotics in the transition regions, near the lines x/>ct = 2 and x/>it = —1/4, 
involves Painleve transcendents. A detailed analysis is presented in [1]. 

2. Main result 

In order to state our main result we first need to recall a few things. 
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Scattering data. Associated with w{x, t) is a self-adjoint Sturm-Liouville operator 

1 / d^ 1\ 

~ w{x,t) V dx2 + 4y ' (2.1) 
D{H{t)) = H'^{R) C L'^{R,wdx). 

Here L2(M, dx) denotes the Hilbert space of square integrable (complex- valued) func- 
tions over M and H'^(U.) the corresponding Sobolev space. 

By our assumption p.5|) the spectrum of H{t) (independent of t) consists of an 
absolutely continuous part [^,oo) plus a finite set of eigenvalues Xj G (0,4^), 1 < 
j < N, possibly empty ( "solitonless case"), see [9l Theorem 2.1]. For j = 1, . . . , we 
denote 

< Kl < K2 < ■ ■ ■ < Kn < - ■ 

Moreover, associated with the continuous spectrum is a (right) reflection coefficient 
R^(k,t), and associated with every eigenvalue Xj is a (right) norming constant 7+j(t) 
(see the next section for details). 

One-soliton solution. Given k S (0, ^) and 7 > we have the corresponding one- 
soliton solution u^J^^"'^ = Uso\{x,t) given by 

,2 



32x/t a(y(x — xct)) 
Uso\{x - xct) = — r— ^-^ , (2.2 



with 



2(i-K^) 
where y{x) is given implicitly by 

l + «(?/)TTi 



= ^e-^'^^ (2.3a) 
^ G(2,cx)), (2.3b) 



^ = y + log . . ,.,'1:2:; • (2-4) 



The momentum w'^^f^'^ = Wso\{x,t) of this solution is given by 

, , 16^2 a(y(x — xct)) \^ , . 

Note that the one-soliton solution has the form of a single peak which is symmetric 
with respect to its center 

1 , 7^ 1 1 + 2k 

^0 = i^^og— + log — . 

2k 2k 1 — 2k 



4 



A. BOUTET DE MONVEL, A. KOSTENKO, D. SHEPELSKY, AND G. TESCHL 



The maximum at xq is given by Usoi^max = ^ iL^k^ ^"^^ taller waves are narrower and 
travel faster (for fixed x). See Fig. [H where Uso\{x — xq) for fixed x and different values 
of K is displayed. 



Figure 1. One-solitons for fixed x = 1 and k = .15, .25, .3 

Note also that if k = k(x) is changing in such a way that the soliton velocity 

2x 



is fixed, then, as x ^ 0, the form of the soliton approaches that of a peakon, fe"'^', 
which is a nonsmooth, weak solution of the Camassa-Holm equation with x = 
(see [4]). See Fig. [21 where nsoi(a; — xq) for fixed velocity v and different values of 
(x, k) is displayed. Note that Usoi.max = v — 2x. 



Uso\ (x) 




Figure 2. Fixed velocity v = 2.5 and (x, k) = (1, .22), (.5, .39), (.1, .48) 

Now we are ready to state our main result. 

Theorem 2.1 (solution asymptotics). Suppose u{x,t) is a classical solution of the CH 
equation (jl.lh satisfying (jl.5p for some integer I > 1. Let 

R{k) = R+{k, 0) and kj, -/j = 7+,i(0), j = 1, . . . , iV, 

be the (right) scattering data associated with H{0) and the initial condition w{x,0). 

There are four sectors in the {x, t) half-plane in which the long-time asymptotics of 
a solution of the CH equation satisfying (jl.Sp are given by the following formulas: 
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(i) The "soliton" region c := > 2 + C for any small C > 0. Let 

^ j = l,...,N, 



' 2(1 -.2)^ 

and let e > sufficiently small such that the intervals [cj — e, cj + e] are disjoint 

and lie inside (2,cxd). 

(ii) If \^ — Cj\ < £ for some j, one has 

u{x,t) = Uj{x - - KCjt) + 0(t"'), (2.6) 
where uj = u^^^^'"'^^ is the one-soliton solution formed from kj and 

N 



i=j+i * 



and having a phase shift 

i=j+i 

For j = N the product is 1 in ()2.7p and the sum is in (j2.8p . 
(12) V \ — Cj\ > £ for all j, one has 

u{x,t) = 0{t"^). (2.9) 
(ii) The "first oscillatory" region 0<c:=-^<2 — C for any C > 0. Here 



I ,s / 2xA:o(c)z^o(c) / 2xfco(c)^ ^ (,\ , , ( \\ 

+ 0(t^") (2.10) 
/or any ^ < a < 1 provided I > 5, where 



I 1 + c-VTT^ 

^o(c) = - Y , (2.11) 

Mc) = log(l - \R{ko{c))\^), (2.12) 
5o{c) = j- arg{R{ko{c))) + avg{T{m{c))) 

. ^ 8xfco(c)^(3/4-fco(c)^) 
-^°^^^'°H (l/4 + ^o(c)^)3 

N 

+ 45]arctan(-^)+- / log(|C - fco(c)|) dIog(l - |i?(C)r 



^•=1 ■A;o(c)^ vr 72(c) 



+ 4.Mc)f:.og(i±|^)+l^M£)/ MffipdC. (2.13, 
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with S(c) = {—ko{c),ko{c)) andT the Gamma function. 



(iv) 

fast decay 



2" oscillatory 



region 



\t 



^ (iii) 
\ 

\ 
\ 
\ 
\ 



(ii) 

1"* oscillatory 
region 




C = Cn 



Figure 3. The different regions of the (x, i)-half-plane, c 



(iii) The "second oscillatory" region —\ + C<c=^<Q for any C > 0, where 



u{x, t) 



2xfco(c)i^o(c) 



sm 



2>cko{cf 



-i-z/o(c) log(t) + 5o(c) 



2xki{c)ui{c) 



sm 

\l{\ + k,{cY){k,{cY-l)t \{\ + k,{cYY 
+ o(t-") 



2xki{cY 



t + vilogit) - 5i{c) 



(2.14) 



LONG-TIME ASYMPTOTICS FOR THE CH EQUATION 

for any h < a < 1 provided I > 5. Here, for £ = 0,1, 



, , , 1 / 1 + c- (-l)Vl + 4c 

kiic) = ^]J , (2.15) 

Mc) = -^logil-\R{ke{c))f), (2.16) 
Se{c) = \- (-l)^arg(i?(fe,(c))) + arg(r(ii.,(c))) 



1 



3=1 

(-!)■ 



[ log(|C-fc,(c)|)dlog(l-|i?(C)|') 

. (A - hr,(A 



fci(c) + /co(c) 

+ 4x.,(c)flogl±^ + l^/ '^^^^^C. (2.17) 

with E(c) = (-oo,-A:i(c)) U (-A;o(c), A;o(c)) U (A;i(c),oo), 0' = 1, araci 1' = 0. 
(iv) The 'fast decay" region c = ^ < — \ — C for any C > 0, where 

u{x,t) = 0{t-^). (2.18) 

We can also give the asymptotics of the momentum w{x,t). 

Theorem 2.2 (momentum asymptotics). Under the same hypotheses and with the 
same notations the momentum of the solution behaves as follows: 

(i) In the soliton region c = ^ > 2 + C . 

(11) If — Cj\ < e for some j, one has 

w{x, t) = Wj{x - - KCjt) + 0(i"') (2.19) 

where wj is the momentum of the one-soliton solution formed from the same 
parameters {Kj,7j} and with the same phase shift as above. 

(12) If \^ — Cj\ > £ for all j, one has 

w{x,t) = >c+0{t-^). (2.20) 

(ii) In the first oscillatory region 0<c=-^<2 — C for any C > 0, one has 



+ 0{t-^), (2.21) 
for any ^ < a < 1 provided I > 5. 
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(iii) In the second oscillatory region —\ + C<c=-^<0 for any C > 0, one has 



3 



+ 0(t-"), (2.22) 

/or any ^ < a < 1 provided I > 5. 
(iv) /n t/ie fast decay region < ~\ ~ C for any C > 0, we have 

w{x,t) = >c+0{t-^). (2.23) 

In particular we recover the fact that a pure soliton solution (i.e., R{k) = 0) asymp- 
totically splits into single solitons with associated phase shifts. This was shown only 
recently by R.S. Johnson |23] (for two solitons) and in the general case by Y. Matsuno 
|29j . For further results on solitons of the CH equation and their stability we refer to 
A. Constantin and W. Strauss [14J, L.-C. Li [26], and K. El Dika and L. Molinet [18] . 

Notice that the oscillatory regions (ii) and (iii) match at x = 0. Indeed, as x — > 
with X < 0, fci — > oo in (j2.15p and thus the amplitude of the second term in (j2.14p 
vanishes, while the parameters of the first term in p.l4p match those in p.lOp . As 
for the transition between the other regions, we have already noticed in Section [1] that 
there exist transition zones, where the asymptotics are described in terms of Painleve 
transcendents. More precisely (details are given in [5), these zones are: 

(trl) \x/Kt - 2| t2/3 < Const, 
(tr2) \x/>tt + 1/4| t2/3 < Const. 

It should also be emphasized that unlike the (modified) Korteweg-de Vries (KdV) 
equation (originally considered in [l6]), the asymptotic form is given implicitly, however, 
to leading order this fact only manifests itself in additional phase shift. More precisely, 
the term in the soliton region (as already pointed out in [29]) and the last two terms 
in (^j(c), have no analog in the (modified) KdV equation (cf. T6], respectively [20]). 
For results on CH on the half-line we refer to A. Boutet de Monvel and D. Shepelsky 

mm- 

Finally, note that if m(x, t) solves the CH equation, then so does u(— x, —t). Therefore 
it suffices to investigate the case t +oo. 



3. The Inverse scattering transform and the Riemann-Hilbert problem 

In this section we derive a vector Riemann-Hilbert problem directly from the scat- 
tering theory for the differential operator ()1.3p . We begin by recalling some required 
results from scattering theory, respectively the inverse scattering transform for the CH 
equation from [HKH] (see also [28]). 
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Recall also that by virtue of the unitary Liouville transform 

f{x)^ f{y)=wix)'/^fix), 

"+~/ lw(r) \ . (3-1) 



y = x- \ 1 dr 



the Sturm-Liouville operator Hit) introduced in (|2.1|) can be mapped to a self-adjoint 
Schrodinger operator 

k\ dy^ 4y (3.2) 

D{H{t)) = H'^{R) C l2(M). 

where 

q(^y t) - ^ ^^^(^'^) _ wix,t) - K _ hKWxix.tf 
' 4 4ii;(x,t) 16 tt;(2;,t)^ 

From our assumption (jl.Sp it follows that q{y,t) G L-'^(M, (1 + |y|)'+-'^dy). 

Lemma 3.1. There exist two Jost solutions ijj±{k,x,t) which solve the differential 
equation 

H{t)i;±{k,x,t) = ^(^^ + k^^ij±{k,x,t), Im(fe) > 0, (3.3) 

and 

lim e=f^'=^V±(fc,2;,i) = 1. (3.4) 

x— >±oo 

Both 7p±{k,x,t) are analytic for Im(A;) > and continuous for lm{k) > 0. For large k 
we have 

1/4 / ^±00 1 / 1 \ \ 

fc(*..,*)=e-(.-^''-.W)^l^(lTy^ ,(,.,t)d.^ + 0(^)) (3.5) 

as k ^ oo, where 



H^^{w) = jj^^^-iyr (3.6) 
is a conserved quantity of the CH equation. 

Proof. This is immediate from the corresponding results for (j3.2p (cf., e.g., [15] or [28]) 
by virtue of our Liouville transform ()3.ip . Just observe 

V^±(A;, X, t) = e-'^^^-i(-) ^ V^±(fc, y, t), 

where 'tp±{k,y,t) are the Jost solutions of (j3.2p . □ 
Furthermore, one has the scattering relations 



r(A;)?/'zF(/c,x,t) = ?/;±(fc,x,t) +/?±(A;,t)V'±(A:,x,t), A; G M, (3.7) 

where T{k), R±{k,t) are the transmission, resp. reflection coefficients. We have sym- 
metry relations 

R±i-k,t) = R±{k,t) and T{-k) = T{k). 
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Note also that if T{k), R^{k,t) are the corresponding quantities for H(t), then T{k) = 
e'*^^-i(^)r(/c), R+{k,t) = R+{k,t), and R-{k,t) = e'^^^^-^^'") R_{k,t) and hence all 
results known for (j3.2p readily apply in our situation. In particular, they have the 
following well-known properties: 

Lemma 3.2. The transmission coefficient T{k) is meromorphic for Im(/c) > with 

simple poles at Iki, . . . , \kn , where as above kj = \J\~ ^^j € (0) f ); ^^'^ ^-^ continuous 
up to the real line. Asymptotically we have 

T{k) = ei^-f^-iW(l + 0{k~^)). (3.8) 
The residues ofT{k) are given by 

ResiK^. T{k) = ifij{t)j+j{t)'^ = ifijjlj, (3.9) 

where 

7_^j(t)-2 = x-i / ij+{iKj,r,tfw{r,t)dr (3.10) 
Jr 

and ip^ {iKj ,x,t)=fij {t)'ip^ {iKj ,x,t). 
Moreover, 



T{k)R+{k,t) +T{k)R_{k,t) = 0, \T{k)\^ + \R±{k,t)\^ = 1. (3.11) 

Note that one reflection coefficient, say R{k,t) = R^(k,t), and one set of norming 
constants, say 7j(t) := 7+j(t), suffices. 
The time dependence is given by (see P] ) : 

Lemma 3.3. The time evolutions of the quantities R{k,t) and ^j{t) are given by, 



R{k,t) = R{k)e , (3.12) 



7i(i)=7,e^^^-'=5 (3.13) 
where R{k) = R{k,0) and 7j = 7j(0). 

Vector Riemann— Hilbert problem. We will set up a vector Riemann-Hilbert prob- 
lem as follows. Let m{k,x,t) = (mi(A;,x,t) m2{k,x,t)^ be defined by 

{^^^y/^ (T{k)ilj^{k,x,t)e'^y il)+{k,x,t)e-''^y^ , Im(A;) > 0, 



{^^Y/^{ij+{-k,x,ty''y T{-k)iP_{-k,x,t)e''^y'^ , lm{k)<0. 



(3.14) 



We are interested in the jump condition of m{k,x,t) on the real A:-axis (oriented from 
negative to positive). To formulate our jump condition we use the following convention: 
when representing functions on M, the lower subscript denotes the non-tangential limit 
from different sides. By mj^{k) we denote the limit from above and by m-{k) the one 
from below. Using the notation above implicitly assumes that these limits exist in the 
sense that m{k) extends to a continuous function on the real axis. In general, for an 
oriented contour S, m^{k) (resp. m-.{k)) will denote the limit of m(K) as k, ^ k from 
the positive (resp. negative) side of S. Here the positive (resp. negative) side is the one 
which lies to the left (resp. right) as one traverses the contour in the direction of the 
orientation. 
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Theorem 3.4 (vector RH-problem) . Let S+{H{0)) = {R{k), (^^,7^), j = 1,...,N} be 

the right scattering data of the operator H{0) associated with the initial data w{x,0). 
Then m{k) = m{k, x, t) defined in (|3.14p is a solution of the following vector Riemann- 
Hilbert problem. Find a function m{k) which satisfies: 

(i) The analyticity condition: 

m{k) is meromorphic away from the real axis with simple poles at ^iKj. 

(ii) The jump condition, for A; G M; 

mj^{k) = m^{k)v{k), 

fl-\Rik)\^ -:R(fc)e-**('=)\ (3.15) 
Vi2(A:)e**('=) 1 



v{k) 

(iii) The pole conditions, for j = 1, . . . , N : 



ReSiK, m{k) = lim m{k) { .^2gt$(iK,) gj • (3-16) 



In (ii) and (iii) the phase is given by 



ji + k^ t 



(iv) The symmetry condition 



(v) The normalization 



m{-k)=m{k)[\ M. (3.18) 



lim m{k) = (1 1). (3.19) 

fe^oo 



Remarks, (a) Note that detf(A;) = 1 and v{—k) = aiv{k) ^ai with fii = (5 J), 
(b) Note also that (iii) and (iv) imply the pole conditions, for j = 1, . . . , N: 

-i72ei*{i«i) 



Hes-^K mik) = lim m(k) , 

Proof. The jump condition (j3.15p is a simple calculation using the scattering relations 
(j3.7p plus p.lip . The pole conditions follow since T{k) is meromorphic for Im(A;) > 
with simple poles at ikj and residues given by (j3.9p . The symmetry condition holds by 
construction and the normalization (jS.lOp is immediate from (|3.5p and 

T{k) = e'*--^-!^'") (^1 + y"^°° q{r, t)dr{2ik)-^ + 0{k-^)^ (3.20) 

which implies 

m{k,x,t) = {l l)+Q+(y,t)-L(l _l)+0(A:-2), (3.21) 

where Q+{y,t) = fy°° q{r,t)dr. □ 
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Observe that the pole condition at i/tj is sufficient since the one at —inj fohows by 
symmetry. Hence, the Riemann-Hilbert problem for the Camassa-Holm equation is, 
for given scattering data 5+, to find a sectionally meromorphic vector function m{k) 
satisfying (j3.15p ~ (j3.19p . We will show that the solution given in the above theorem 
is in fact the only one in Corollary 13.91 below. Moreover, it should be pointed out 
that except for the phase, this Riemann-Hilbert problem is identical to the one for the 
Korteweg-de Vries equation (cf. [201 Thm. 2.3]). 

Next we note the following useful asymptotics 

Lemma 3.5. The function m{k, x,t) defined in (|3.14p satisfies 

^I^lMA = e-^ (3.22) 
m2{^,x,t) 

and 

mi{k, X, t)m2{k, X, t) = ^355(^1 + ^u{x, t) {k - ^) + 0{k - . (3.23) 

Proof. The Jost solutions admit the representation 'ip±{k, x, t) = e^^^^g±{k, x, t), where 
g±{k,x,t) are the solutions of the integral equations 

U2 I 1 r±oo 

g^{k,x,t) = l±^-j^ J {l-e'f^'''^'^-'-^)g±{k,r,t)iwir,t)->,)dr. (3.24) 

Since w satisfies ()1.5p . the solution of (j3.24p exists and is unique (see, for example, 
[28]). Moreover, g± is analytic for lm(k) > 0. 
Since fe^ + 1 = (A; - i)(/c + |), we get 

g±{k,x,t) = l±-{k- ^)F±{x,t) + 0{k- ^f, k ^ i/2, (3.25) 

where 

F±ix,t) = / (e±(^-^) - l){w{r,t) - x)dr. (3.26) 

J X 

Moreover, differentiating with respect to x, we see 

g'±ik,x,t) = ±-ik - |)F4(x,t) + 0{k - i)2, (3.27) 

with 

/"±oo 

F4(x,i) = ± / e±(^-'')(u;(r,f) - x)dr. (3.28) 

J X 



Using 



u{x,t) = (1 - a^) ^ (w{x,t) - x) = - I Q-\^-''\(w(r,t) - x)dr, (3.29) 

2 ./H 



we thus obtain 



i a2 



V'+(A;,a;,t)V'-(A:,x,i) = 1 + -(F+(a;, t) - F_(x, t))(A; - |) + 0(A: - ^) 

0^ 



1 + -{2u(x, t) - Hoiu)){k - I) + 0(A: - ^f, (3.30) 

0^ 
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where 

Hq{u) = / {w{x,t) — M:)dx = / u{x,t)dx (3.31) 
Jm Jr 

is a conserved quantity of the CH equation. 

Furthermore, straightforward calculations show that 



T{k) 



2ik 
i 



1 + - (F+(x, t) - F_ (x, t) - F'^{x, t) - F'_ (x, t)) {k-^) + 0{k - \f 

i \2 



l--HQ{u){k-'^) + 0{k-^)\ (3.32) 



where W(/, g) = fg' — f'g is the usual Wronskian. Therefore, 

T{k) = 1 + -Ho{u) {k-'^)+ 0{k - i)^ . (3.33) 



Substituting (l3:25]l and (13331) into (IXTiD . we arrive at IK22\\ . Substituting (13301) and 
(I333D into we obtain iK23\f . □ 



Regular Riemann— Hilbert problem. For our further analysis it will be convenient 
to rewrite the pole condition as a jump condition and hence turn our meromorphic 
Riemann-Hilbert problem into a holomorphic Riemann-Hilbert problem. 

Choose e so small that the discs \k — iKj\ < e lie inside the upper half plane and do 
not intersect. Then redefine m{k) in a neighborhood of iKj, resp. —inj according to 

I k ^ J I ^ ? 

(3.34) 

|A; + iKjl < e, 
else. 

Note that we redefined m{k) such that it respects our symmetry ()3.18p . 

Then a straightforward calculation using Resi^ m{k) = limfc_^iK(A; — m)m{k) shows: 

Lemma 3.6 (regular RH-problem). Let Cj be the circle \k — iKj\ = e with e > as 
above, I < j < N . Let m[k) be defined as in (|3.34p . Then m{k) is a solution of the 
following vector Riemann-Hilbert problem. Find a function rh{k) which satisfies: 

(i) rh[k) is holomorphic away from the real axis and from the circles Cj and Cj, for 
j = l,...,N. 

(ii) The jump condition p.lSp across the real axis. 



m(k) 



rh[k) 



m{k) I fc+iKj 



,0 



1 



m{k), 



14 A. BOUTET DE MONVEL, A. KOSTENKO, D. SHEPELSKY, AND G. TESCHL 

(iii) The additional jump conditions across the circles Cj, Cj, for j = 1, . . . ,N: 

( 1 o\ 

rh+{k) = m^{k) i^V^^^'^j' , k e Cj, 

\ k—iKj ^/ 

k+iKj , k ^ Cj, 
1 J 

where Cj is oriented counterclockwise and Cj is oriented clockwise. 

(iv) The symmetry condition ()3.18p . 

(v) The normalization condition ()3.19p . 



(3.35) 



6^^ 



Figure 4. Contour of the regular RH problem 



Uniqueness result. Next we turn to uniqueness of the solution of this vector Riemann- 
Hilbert problem. This will also explain the reason for our symmetry condition. We 
begin by observing that if there is a point ki G C, such that m(fci) = (O O), then 
n{k) = ■^^^m{k) is a solution of the associated vanishing Riemann-Hilbert problem, 
i.e., it satisfies the same jump and pole conditions as m{k) but the normalization 
now reads lim^^oo m{iK) = (O O) . In particular, there is a whole family of solutions 
m{k)+'dn{k) for any G C. However, these solutions will clearly violate the symmetry 
condition unless = 0! Hence, without the symmetry condition, the solution of our 
vector Riemann-Hilbert problem will not be unique in such a situation. Moreover, a 
look at the one-soliton solution verifies that this case indeed can happen. 

Lemma 3.7 (one-soliton solution). Suppose that the reflection coefficient vanishes, i.e., 
R{k,t) = and that there is only one eigenvalue k £ (0, ^), with the norming constant 
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■j{t). Then the unique solution of the Riemann-Hilbert problem p.l5p - (j3.19p is given 
by 

moik) = {f{k) fi-k)) (3.36) 
fik) = -^(l + a'^'^ 



In particular, 



1 + a V k — IK 



2k 



^) = (1 _ 4^2)2 "(^' ^) \^(^ + ^)) + 1 _ 4^2 "(^'^) ) ' (3-37) 
w;(x, t) = X 1 + , (3.38) 



1 -4k2 (1 + a(y,t))2 

i + a(y,t)S 

Proof. By assumption the reflection coefficient vanishes and so the jump along the real 
axis disappears. Therefore and by the symmetry condition, we know that the solution 
is of the form mo{k) = {f{k) f{—k)) where f{k) is meromorphic. Furthermore the 
function f[k) has only a simple pole at Ik, so that we can make the ansatz f{k) = C + 
I?|±l^. Then the constants C and D are uniquely determined by the pole conditions and 
the normalization. Formulae (|3.37p - p.38p are obtained applying Lemma 13.51 formula 
(13:23]) . □ 

In fact, observe that f{ki) = f{—ki) = if and only if /ci = and 2k = ^^gt^CiK)^ 
Furthermore, even in the general case m{ki) = (O O) can only occur at /ci = as the 
following lemma shows. 

Lemma 3.8. If m{ki) = (O O) for m defined as in (I3.14|) . then ki = 0. Moreover, 
the zero of at least one component is simple in this case. 

Proof. By (I3.14P the condition m{ki) = (O O) implies that the Jost solutions 'ip-{k,x) 
and 'ilj+{k,x) are linearly dependent or that the transmission coefficient T{ki) = 0. 
This can only happen at the band edge, fei = or at an eigenvalue ki = ikj. 

We begin with the case ki = ikj. In this case the ip-{k,x) and ipj^{k,x) are linearly 
dependent. Moreover, T{-) has a simple pole aX, k = ki since the derivative of the 
Wronskian W(/c) = ipj^{k,x)^'_{k,x) — 'ip'_^_{k, x)'ip-{k, x) does not vanish by the well- 
known formula 

^\N{k)\k=kr = -— f Tp+{ki,x)ip.{ki,x)w{x)dx ^ 
dk X 

(cf. Lemma l3.2p . The diagonal Green's function g{z,x) = \N{k)'^^ip^{k, x)ip^{k,x) is 
Herglotz as a function of z = — A;^ and hence can have at most a simple zero cit z — — /c^. 
Since z ^ —k^ is conformal away from z = the same is true as a function of k. Hence, 
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if 'ip+{iKj,x) = ^-{iKj,x) = 0, both can have at most a simple zero at A; = ikj. But 
T{k) has a simple pole at iKj and hence T{k)^-{k, x) cannot vanish at A; = ikj, a 
contradiction. 

It remains to show that one zero is simple in the case ki = 0. In fact, one can show 
that 

Aw(A;)U=fc, /O 

in this case as follows: first of all note that ip±{k) (where the dot denotes the derivative 
with respect to k) again solves 

Hi;±{ki) = ^(^^ + kl^ V'±(A:i) 

if ki = 0. Moreover, by \N(ki) = we have ■i/'+(^i) = ctj)-{ki) for some constant c 
(independent of x). Thus we can compute 

W(fci ) = W(V3+ (fci ) , V^- (fci ) ) + W(V'+ ( A;i ) , ^_ ( A;i ) ) 

= c-i W(V'+(A:i),V+(A:i)) +cW(V'-(A:i),V'-(A:i)) 

by letting x +oo for the first and x — oo for the second Wronskian (in which case 
we can replace ip±{k) by e^'*^^), which gives 

W(A;i) = -i(c + c-i). 

Hence the Wronskian has a simple zero. But if both functions had more than simple 
zeros, so would the Wronskian, a contradiction. □ 

By [20i Theorem 3.2] we obtain 

Corollary 3.9. The function m{k, x, t) defined in (I3.14p is the only solution of the 
vector Riemann-Hilbert problem (I3.15p - (l3.19p . 

4. Conjugation and Deformation 

This section demonstrates how to conjugate our Riemann-Hilbert problem (with 
respect to the augmented contour) and how to deform our jump contour, such that 
the jumps will be exponentially close to the identity away from the stationary phase 
points. Throughout this and the following section, we will assume that the R{k) has 
an analytic extension to a small neighborhood of the real axis. This is for example the 
case if we assume that our solution is exponentially decaying. This assumption can 
then be removed using analytic approximation. 

For easy reference we note the following result: 

Lemma 4.1 (conjugation). Let T, be a part of some contour S. Let D be a matrix of 
the form 

where d: C\S— >C is a sectionally analytic function. Set 

m{k) = m{k)D{k), (4.2) 



LONG-TIME ASYMPTOTICS FOR THE CH EQUATION 



17 



then the jump matrix transforms according to 



v{k) = D^{k)-^v{k)D+{k). 



(4.3) 



If d satisfies d{—k) = d{k) 



and lim, 



■A;— »oo 



d{k) = 1, then the transformation 



■m{k) = m{k)D{k) 



respects our symmetry condition, that is, rh{k) satisfies (j3.18p if and only ifm{k) does, 
and our normalization condition. 

In particular, we obtain 



In order to analyse the regular vector RH problem from Lemma [3.61 there are two cases 
to distinguish. 

(a) If ^{iKj) < then the corresponding jump matrix (j3.35p is exponentially close to 
the identity as t — > +00 and there is nothing to do. 

(b) Otherwise we use conjugation to turn the jumps into one with exponentially de- 
caying off-diagonal entries. 

It turns out that we will have to handle the jumps across Cj and Cj in one step in 
order to preserve symmetry and in order to not add additional singularities elsewhere. 

Lemma 4.2. Assume that the Riemann-Hilhert problem for m has jump conditions 
near in and —ik given by 



Then this Riemann-Hilbert problem is equivalent to a Riemann-Hilbert problem for 
rh = mD which has jump conditions near in and —ik given by 




(4.4) 




(4.5) 




(4.6) 
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and all remaining data conjugated by 



D{k) 






k+\i 
k—\n, 





m\ < e, 



\k + ikI < £, 



(4.7) 



else. 



The jump along the real axis is of oscillatory type and our aim is to apply a contour 
deformation following [16] such that all jumps will be moved into regions where the 
oscillatory terms will decay exponentially. Since the jump matrix v contains both 
exp(t<I>) and exp(— 1$) we need to separate them in order to be able to move them to 
different regions of the complex plane. 

We recall that the phase of the associated Riemann-Hilbert problem is given by 

^k) = -iY^^ + 2ik-. (4.8) 



t 



Let 



c{y,t) = — . 

xt 



(4.9) 



The stationary phase points, i.e., ^'{k) = 0, are given by it/co and =b/ci, where 



ko = koic) 



1 + c- Vl + 4c 



1 / 1 + c + a/1 + 4c 



(4.10) 



ki = ki{c) = - ^ 

There are four cases to distinguish: 
(i) 2x < J, i.e. c > 2. In this case 
ko, ki M, 

> Re(<^) < for Im(A;) > near R U i(0, kq) and 

> Re(<l>) > for Im(A;) > near i(Ko, ^), where 



i./l-i 



(4.11) 



for J < 2>c for notational convenience later on. 
< f < 2x, i.e. < c < 2. 



(iii) 



We will set kq 

In this case 

> ko £R, ki^ R, 

> Re(<&) > for Im(A;) > near (-A:o, ko) U i(0, i) and 
Re(<I>) < for lm{k) > near (-oo, -ko) U (ko, oo). 
-J < I < 0, i.e. -1/4 < c < 0. In this case 

l> ko, ki G R, 

t> Re(^>) > for Im(A;) > near (-oo,-A;i) U {-ko,ko) U (/ci, oo) U i(0, i) 

> Re(<I>) < for Im(A;) > near {-ki, -ko) U {ko, ki). 



and 
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(iv) J < — J, i.e. c < — 1/4. In this case 

> ko, fci M and 

> Re($) > for Im(/c) > near MU i(0, i). 
The situation is depicted in Figure [5l 



c < -1/4 



-1/4 < c < 



< c < 2 



2 < c 

















5 








+ 





0, 



0" 



+ 



Figure 5. Sign of Re($(A;)) for different values of c 



>c t 



Accordingly we will introduce 



c < -4, 



(-00, -/ci) U {-ko, ko) U (A;i, 00), < c < 0, 

(-A:o,A;o), < c < 2, 

0, 2 < c. 



(4.12) 



As mentioned above we will need the following factorizations of the jump condition 

v{k) = b^{k)-^b+{k), 

where 

'1 i?(fc)e-**('=)' 



b^{k) = 
for k G M \ S(c) and 







1 



6+ (A;) 



v(A;) = B_{k) 



1 n-\R(.k)y 



i-i/?(fc)i 



1 



where 



i?_(A;) 



1 0^ 

' l-|/?(fc)P y 



R{k)e 



-t*{fc) ' 



,0 



l-l/?(fc)l^ 



1 



(4.13) 
(4.14) 



(4.15) 



(4.16) 



for k G S(c). 

To get rid of the diagonal part in the factorization corresponding to A; E 5](c) and to 
conjugate the jumps near the eigenvalues we need the partial transmission coefficient 
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with respect to c defined by 



T{k,c) 



n 

K0<Kj<l/2 

^ k+iK 



C> 2, 



n k-iKi ( 2ni It, 



log(inc)P) 

(c) c-fc 



dC), c<2, 



(4.17) 



for /c G C \ S(c). Thus T{k,c) is meromorphic for k £ C \ S(c). Note that T{k,c) can 
be computed in terms of the scattering data since |r(A;)p = 1 — \R^{k,t)\'^ . Moreover, 
we set 



ri(c) = iog(r(i,c)) = < 



c> 2, 



K0<Kj<l/2 

f log + i r dr c < 2 



(4.18) 



Note that combining r(A;) = e'^^-^^'"'>T{k,c) for c < with (f^IMj) shows 

^ 1 + 2K,- , 2 /"logdncp^^^ 



i7_,(.)=2j:iogp|^ + ^/^ 



+ 4C^ 



(4.19) 



Theorem 4.3. The partial transmission coefficient T{k,c) is meromorphic in C\S(c) 
with simple poles at iKj and simple zeros at —iKj for all j with kq < Kj, and satisfies 
the jump condition 

T+{k,c) =T_{k,c){l-\R{k)\'^), forke^ic). (4.20) 

Moreover: 

(i) T{-k, c) = T{k,c) -\ k£C\ S(c). 

(ii) T{—k, c) = T{k, c), /c G C, in particular T(k, c) is real for k G iM. 

(iii) If c < 2 the behaviour near k = is given by T{k,c) = T{k){C + o(l)) with C 7^ 
for Im(/c) > 0. 

Proof. That mj are simple poles and — i^j are simple zeros is obvious from the Blaschke 
factors and that T(A;, c) has the given jump follows from Plemelj's formulas. Properties 
(i), (ii), and (iii) are straightforward to check. □ 

Now we are ready to perform our conjugation step. Introduce 



D{k) 



1 


■ 2 tt?(lKA 
"JjC ^ 3 ' 



17, e ^ ^' 



Do{k), \k — IKjl < e, Kq < Kj, 



k+i 



'^^ \ Do{k), \k + IKjl < e, Kq < Kj, 



[Do{k), 



else, 



where 



Do{k) 



T{k,c)-^ 
T{k,c) 
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Observe that D{k) respects our symmetry: 

Now we conjugate our problem using D{k) and set 

rh{k) = m{k)D{k). 
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(4.21) 



Note that even though D{k) might be singular at A; = (if c < 2 and -R(O) = — 1), fn{k) 
is nonsingular since the possible singular behaviour of T{k,c)^^ from Do{k) cancels 
with T{k) in m(fc) by virtue of Theorem 14.31 (iii). 

Then using Lemma l4.ll with S = S(c) and Lemma 14.21 the jump corresponding to 
kq < Kj (if any) is given by 



v{k) 



v{k) 



i^2(,t*(iKi)j.(fc^c)-2 



I J\i \Ki j I — S ^ 



1 



^ i72e«*(''=j)T(fc,c)2 

and corresponding to kq > Kj (if any) by 



(4.22) 



\k + IK 



'J I 



v{k) 



y2e**('''j'T(fc,c)-2 

iT,2e**(-j)T(fc,c)2' 
1 



\k — IK 



'j I ~ ^' 



\k + iKjl = e. 



(4.23) 



In particular, all jumps corresponding to poles, except for possibly one if kj = kq, are 
exponentially close to the identity. In the latter case we will keep the pole condition 
for Kj = Kq which now reads 



ReSiK, rh^k) = lim m{k) • 2„t$(iK 



k—fin 











ij^em^j)T{[Kj,cy 



Res.,..rnik)= lim ^k) ('^ -i7|e-^-^^n-„ c)-y 



(4.24) 



Furthermore, the jump along M is given by 

v{k) 

where 



b^{k)-^b+{k), k^i:{c), 
B-{k)-^B+{k), k G S(c), 







1 



1 

iJ(fc)e**W 



T'(fc,c)2 



1 



(4.25) 



(4.26) 
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and 



B-ik) 



,0 



r_(fc,c)- 

' l-\R{k) 

l-\Rik)\ 



T_(fc,c) 



r+(fc,c) p. 
'T+{-fc,c)-^'^ 



-A:)e 



-t$(fc) 



Here we have used 



k G S(c), 



i?(-A;) = R{k), 
T±{-k,c)=T^( 

and the jump condition (j4.20p for the partial transmission coefficient T{k, c) along S(c) 
in the last step. This also shows that the matrix entries are bounded for k G M. near 
k = since T±{-k,c) = T±{k,c). 

Since we have assumed that R{k) has an analytic continuation to a neighborhood of 
the real axis, we can now deform the jump along R to move the oscillatory terms into 
regions where they are decaying. 

There are four cases to distinguish: 

Case (i): c > 2. We set Sj- = {k £ C \ Im(A;) = ite} for some small e such that Sj- 
lies in the region with it Re(<l>(A;)) < and such that the circles Cj, Cj around Hkj lie 
outside the region in between S_ and E+. Then we can split our jump by redefining 
fh{k) according to 

m{k)b+{ky^ , < Im(A;) < e, 

m{k)b^{k)-^, -e<lm{k)<0, (4.27) 
m{k), else. 

Thus the jump along the real axis disappears and the jump along E± is given by 

'h+{k), k£j:+ 

All other jumps are unchanged. By construction the jump along T,± is exponentially 
close to the identity as t — > oo. 

Cases (ii) and {in): < c < 2, respectively — 1/4 < c < 0. We set S± = U 
according to Figure [6] respectively Figure [7] again such that the circles around lie 
outside the region in between S_ and Again note that respectively lie in 
the region with ibRe($(A;)) < 0. 

Then we can split our jump by redefining m{k) according to 

and 

and EL, 



m{k) 



v{k) 



(4.28) 



m{k) 



m{k)b+{k)-^, 
m{k)b-{k)-^, 
m{k)B+{k)-\ 
m{k)B^(k)-^, 
m{k), 



k between 
k between 
k between 
k between 
else. 



and E^, 
and E?_, 



(4.29) 



LONG-TIME ASYMPTOTICS FOR THE CH EQUATION 23 

One checks that the jump along R disappears and the jump along T,± is given by 













A; G SL 















m = {^\' ' , (4-30) 



All other jumps are unchanged. Again the resulting Riemann-Hilbert problem still 
satisfies our symmetry condition (j3.18p and the jump along S± \ {±ko,±ki} is expo- 
nentially decreasing as t — > cxd. 

Case (iv); c < —1/4. We set S± = {A; G C | Im(A;) = ite} for some small e such that 
T,± lies in the region with ibRe($(A;)) > and such that the circles around ±1^-,- lie 
outside the region in between S_ and Then we can split our jump by redefining 
rh(k) according to 

^m{k)B+{k)-'^, 0<lm{k)<e, 

m{k)B^{k)~\ -e<Im(A;)<0, (4.31) 
m{k), else. 
Thus the jump along the real axis disappears and the jump along Ej- is given by 



m{k) 



B+{k), A; G Sh 
B-{k)-\ ke^- 



All other jumps are unchanged. By construction the jump along Ej- is exponentially 
close to the identity as f — > oo. 

Note that in all cases the resulting Riemann-Hilbert problem still satisfies our sym- 
metry condition ()3.18p . since we have 

b±{-k)=(^^ J), ^±(-A;)=(j J)- (4.33) 

In Cases (i) and (iv) we can immediately apply Theorem lA.ll to in as follows: 



Re$>0 


Re*<0 




Rc'I'>0 


Re*<0 


Re*>0 
\ ^\ 




/ ' \ 


/-feo 










V . /K y\ . y 






si 










Re$<0 


Re<I>>0 




Rc4><0 


Re $>0 


Rc$<0 



Figure 6. Deformed contour for— l/4<c<0 
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Re4><0 



Si 



/'Re<I>>0 \ 



Re*<0 

Si 



-fco 



Re4>>0 



v,Re$<0, 



Si 
Re'I«>0 



Figure 7. Deformed contour for < c < 2 

Proof of Theorem I2.1H2.2I (iv). Since v{k) = I + 0(t^') for any I as t ^ oo, the same 
is true for fh{k) = (l l) + ') by Theorem lA.ll (for the case 7 = 0). Hence 

mi{k)m2{k) = mi{k)m2{k) = 1 + 0(t~'), 

= e2Ti(c)^l(i) ^ ^2T,(c) + o(^-^)^ (4 34) 

for k near ^ and the claim fohows from Lemma [3.51 in case the reflection coefficient has 
an analytic extensions. 

Otherwise one has to split the reflection coefficient into an analytic part plus a small 
remainder. One can literally follow the argument of |20l Lemma 6.1]. □ 

Proof of Theorem I2.1H2.2I (i) . If | j — Cj | > e for all j we can choose 7 = in Theo- 
rem [X?TJ Hence as in the proof of (iv), 

mi{k)m2{k) = fhi{k)m2{k) = 1 + 0(t^'), 



m2{^) m2{^) 
Hows as before. 

Otherwise, if || — Cj\ < e for some j, we choose 7 = jj{x,t). Again we conclude 



for k near ^ and the claim follows as before 



mi{k)m2ik) = mi{k)fh2{k) = f{k)f{-k) + 0{t-'), 

!^ ^ g2T,(c)^i(|) ^ ^2T,ic)Ml + o(t-'), (4.36) 
"^2(^) m2(^) /(-^) 

where f{k) is the one-soliton solution from Lemma 13.71 □ 

In the cases (ii) and (iii) the jump will not decay on the small crosses containing the 
stationary phase points and we need to continue the investigation of this problem in 
the next section. 
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5. Reduction to a Riemann-Hilbert problem on a small cross 

In the previous section we have seen that for — 1/4 < c< 2 we can reduce everything 
to a Riemann-Hilbert problem for fh{k) such that the jumps are exponentiahy close 
to the identity except in small neighborhoods of the stationary phase points ztfeg and 
ib/ci. Hence we need to continue our investigation of this case in this section. 

Denote by S'^(ib/c£), i = 0,1 the parts of S+US_ inside a small neighborhood of ±ki. 
We will now show how to solve the two problems on the small crosses ^^^{ki) respectively 
T,^{—ki) by reducing them to Theorem IA.3I This will lead us to the solution of our 
original problem by virtue of Theorem IA.2[ 

Now let us turn to the solution of the problem on 

S"(A;^) = (S+ U S_) n {A: I |A; - ke\ < e} 

for some small e > 0. We can also deform our contour slightly such that T,'^{k£) consists 
of two straight lines. Next, abbreviate 



(l/4 + A:2)2' 

- (-1) —— - (-1) (i/4 + fc|)3 ' (^-1) 

where > for -1/4 < c < 2 and > for -1/4 < c < 0. 
As a first step we make a change of coordinates 

Q = {-lY^;{k-ke), k = ke+^-^Q (5.2) 

such that the phase reads <^>{k) = (-1)^ i ($^ + + 0(C^)). 

Next we need the behavior of our jump matrix near kg, that is, the behavior of T{k, c) 
near k£. 

Lemma 5.1. We have 



1^ {-^] T(k,c), -l/4<c<0, 
XfeS) ^(^'^)' 0<c<2, 

where f£ = —^log{\T{k£)\) > and the branch cut of the logarithm is chosen along the 
negative real axis. Here 

f(fc,c) = n|^exp(--L / log(|fc-C|)dlog(|T(C)P)). (5.4) 
j-J^ k - iKj V 27ri Js(c) / 

The function T{ - ,c) is Holder continuous of any exponent less than 1 at the stationary 
phase points k = kg and satisfies \T{ki,c)\ = 1. 

Proof. This is a straightforward calculation. Holder continuity of any exponent less 
than 1 is well-known (cf. [30]). □ 
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If k{Q) is defined as in (j5.2p and < a < 1, then there is an L > such that 



TikiCi),c) 



<mr, 



where the branch cut of Q ^ is chosen along the negative real axis and 



-11 



k-kf\-i-iYm 



We also have 
Set 



k + ke 

\R{k{Ci)) - R{ke)\ < L\C 



Tik,c) 



k=ki 



(5.5) 



(5.6) 



(5.7) 



r, = i?(A;,)f,(c)-2e(-i)'''^^'°s(4'^'*") 



(5.8) 

and note ue = — 2^1og(l — \R{kf)\'^) since |r^| = \R{k£)\. Then the assumptions of 
Theorem IA.3I are satisfied with r = tq near for S^(/co)- Similarly, for 'E,^{ki) they are 
satisfied with r = n after a conjugation with ( 5 o )• 

Therefore we can conclude that the solution on S'^(/c^) is given by 



M^{k) = I + 
= 1 + 



1 i /_0_ 
1 



-I3i 



^(A; - ke) \Pi 



+ o(t-": 

0_ -(3e 




+ o{t- 



1/2 <a< 1, 



(5.9) 



where is given by 



Pi 



;/^e^('r/4-arg(ro)+arg{r(ii/o))-t$o)^-ii/0 

^^gi(-7r/4-arg(ri)-arg(r(ii/i))+t$i)^ii/l 

i{-7r/4-arg(i?,(A:i))+2arg(fi(c))+i^i log{Akf<!>'{)-aT:g{T{iui))+t^i)^ivi^ 



(5.10) 



(5.11) 



We also need the solution M^{k) on T,^{—k£). We make the following ansatz, which is 
inspired by the symmetry condition for the vector Riemann-Hilbert problem, outside 
the two small crosses: 



Mm 



1 

1 



{-lY 



1 

1 



pe 

^{k + ke) \-Pe 
Now we are ready to finish the proof of Theorem I2.m2.2i 



+ 0(r 



(5.12) 



Proof of Theorem 12.1112.21 (ii) . By Theorem IA.2I we infer 

1 i „ , 1 



m{k)={l 1) + 

+ o(t-") 



k-ko 



iPo -Po) 



/^k + ko 



-Po Po) 



(5.13) 
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and thus 



m(i) = (l 1) + 
8i 



Re 



/3o 



Re 



/3o 

(l-2iA:o)^ 



-2ifco l+2ifco 

ft 

(l+2iA:o)^ 



+ 0(t-") 



+ 0(r 



(5.14) 



Hence 



mi(k)m2{k) =mi{k)fh2{k) = 1 — ■ 
/ 1 4A:oi 



1/4 + kl 



Im(/3o) + 0(r 



/¥[Jt(l/4 + A;2)2 



Im(/3o) + 0(t-") (A: - i) + . , 



(5.15) 



for k near ^ and Lemma 13.51 implies 

w{x^ t) = K ■ 

u{x, t) = — - 



8A:o 



1/4 + A:2 

Tip (174+^ 



Im(/3o) + 0(r 
Im(/3o) + 0{t- 



(5.16) 



Similarly, 

mi(i) 



-,2T] 



m2( 



^2(5) 



.2ri(c) 



1 + 



^ 1/4 + A:2 



Re(/3o) +0(r"), 



implies 



x-y = 2Tiic) + 



Re(/3o) + 0(r 



(5.17) 



1/4 + A:2 

Thus the claim follows in case the reflection coefficient has an analytic extensions. 
Otherwise one has to split the reflection coefficient into an analytic part plus a small 
remainder using [201 Lem. 6.2 and 6.3]. Again one can literally follow the argument 
given there. □ 

Proof of Theorem I2.1H2.2I (iii) . This follows as in the previous case using 
1 i . N 1 



1 i 



Po /9o) 



+ o(t-"). 



(5.18) 



and hence 

w{x^ t) = X ■ 

u{x, t) = — - 



8ko 



— : TT Im(/3n) 

1/4 + A:2 ^'^"^ 



8k^ 



, ^Im(/3i) + 0(t-") 

$'/tl/4 + A;2 ' 



2kn 



^t{i/^ + klY 



Im(/3o) 



2k^ 



^t{l/^ + klY 



Im(/3i) + 0(t-") (5.19) 
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respectively 

□ 

Appendix A. Some results for Riemann-Hilbert problems 

Irr this section we state the results required to prove or main theorems. We will 
assume that S is a nice contour, say a finite number of smooth oriented finite curves in 
C, which intersect at most finitely many times with all intersections being transversal. 
Moreover, suppose the distance between S and {iy \ y > yo} is positive for some yo > 0. 

A.l. Riemann-Hilbert problem for the soliton region. The first result is needed 
in the soliton region. Consider the Riemann-Hilbert problem of finding a function m{k) 
satisfying 

(i) m{k) is sectionally meromorphic with simple poles at zti/t ^ E, 

(ii) the jump condition 

m^{k) = m-{k)v{k), k&T,, 

(iii) the pole condition 

ReSiKm(A;) = lim m{k) ( n 



(iv) the symmetry condition 

m{—k) = m{k) 

(v) the normalization condition 



1 

1 



lim m{k) = (l l) • 

fe^oo 

Clearly the symmetry condition implies that our jump data (E, v) should be symmetric 
as well, that is, 

<-fc)=(? J)' ^eE. (A.l) 

and E is invariant under k ^ —k and oriented such that under the mapping k ^ —k 
sequences converging from the positive sided to E are mapped to sequences converging 
to the negative side. Moreover (ii) and (iii) imply 



/O -i72\ 
Kes-iKm{k) = lim m{k) „ „ 

k-^—'iK yU U J 



Finally we will assume — I||2 < oo and \\v — I||oo < oo. 

We are interested in comparing this Riemann-Hilbert problem with the one-soliton 
problem (where v = 0) in the case where Hf — I||oo and ||f — ]I||2 are small. For such a 
situation we have the following result: 
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Theorem A.l ([20l[2l]). Assume v = v{t) satisfies 

\\vit)-I\\^<pit), 

\\vit)-Ih<p{t) ^ ■ ' 

for some function p{t) ^0 as t ^ oo. Then the above Riemann-Hilhert problem has 
a unique solution for sufficiently large t and the solution differs from the one-soliton 
solution by 0{p{t)) uniformly in k away from S U {±i/t}. 

A. 2. Riemann Hilbert problem for oscillatory regions. For the oscillatory re- 
gions we will need the following result which allows us to reduce everything to a model 
problem whose solution will be given below. 

Theorem A. 2 ( [20\ I24j). Consider the vector Riemann-Hilbert problem 

mj^{k) = m^{k)v{k), G E, 

lim m(k) = (1 1) , (^-^^ 

fe^oo 

with det(f) / and let < a < (3 < 2a, p{t) — > oo 5e given. 

Suppose that for every sufficiently small e > both the and the L°° norms of 
V are 0{t~^) away from some e neighborhoods of some points fcj G S, 1 < i < n. 
Moreover, suppose that the solution of the matrix problem with jump v{k) restricted to 
the £ neighborhood of ki has a solution which satisfies 

M,{k) = I + + 0{p{tr^), \k-k,\>e. (A.4) 
Then the solution m{k) is given by 

1 " M- 

m{k) = {l + l)E^ + 0(/'(r")' (A-5) 

where the error term depends on the distance of k to S . 

Model problem on a cross. Introduce the cross S = Si U • • • U S4 (see Figure [5]) by 
Si = {ue-'^/^ I n > 0}, S2 = {ne'^/^ | n > 0}, 
S3 = {ue^i-/^ \u>0}, S4 = {ue-=^^"/^ I n > 0}. ^^'^^ 

Orient S such that the real part of k increases in the positive direction. Denote by 
B = {(^ G C I Id < 1} the open unit disc. Throughout this section (^^'^ will denote the 
function 6^'^'°^^''^ where the branch cut of the logarithm is chosen along the negative 
real axis (— cxd, 0). 

Introduce the following jump matrices {vj for C G Sj, j = 1, . . . , 4) 

_ A -i?i(c)C''"e-**(?)\ _ / 1 0^ 
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and consider the RHP given by 

M+(C) = Af_(CH-(C), CeS,-, j = 1,2,3,4, 
M(C) ^ I, C ^ oo. 



(Ai 



The solution is given in the fohowing theorem of P. Deift and X. Zhou [16', Sect. 3-4] 
(for a proof of the version stated below see H. Kriiger and G. Teschl [25, Theorem A.l]). 

Theorem A. 3 ([E]). Assume there is some po > such that Vj{() = I for \(\ > po 
and j = 1, . . . ,4. Moreover, suppose that within |C| < po the following estimates hold: 

(i) The phase satisfies $(0) = i^>o G iM, $'(0) = 0, $"(0) = i, and 



±Re($(C))>-|Cl', with 



+ /orCGSiUSg, 
— else, 



HO - m 



<c\cf. 



(A. 



(A.IO) 



(ii) There is some r E D and constants {a, L) G (0, 1] x (0, oo) such that Ri, . . . , R4 
satisfy Holder conditions of the form 



\RiiO-r\<L\C\' 



\R2{0-r\<L\Cr, 



1 



<L\cr, 



<L\cr. 



(A.ll) 



Then the solution of the RHP (jA.Sh satisfies 

1 i /O 



for I CI > pq, where 



/3 



'I'e 



i{7r/4-arg(r)+arg(r(i!^)))g-it*o^-ii' 



-— log(l - |r|2). 
zvr 



(A.12) 



(A.13) 
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Furthermore, if Rj{C,) and $(C) depend on some parameter, the error term is uniform 
with respect to this parameter as long as r remains within a compact subset of B and 
the constants in the above estimates can be chosen independent of the parameters. 
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